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ABSTRACT 

We analyse the B — Rc colors of galaxies as functions of luminosity and local galaxy density using a 
large photometric redshift catalog based on the Red-Sequence Cluster Survey. We select two samples of 
galaxies with a magnitude limit of Mn^ < —18.5 and redshift ranges of 0.2 < z < 0.4 and 0.4 < z < 0.6 
containing ^ 10^ galaxies each. We model the color distributions of subsamples of galaxies and derive 
the red galaxy fraction and peak colors of red and blue galaxies as functions of galaxy luminosity 
and environment. The evolution of these relationships over the redshift range of z^O.5 to Z'^0.05 is 
analysed in combination with published results from the Sloan Digital Sky Survey. We find that there 
is a strong evolution in the restframe peak color of bright blue galaxies in that they become redder 
with decreasing redshift, while the colors of faint blue galaxies remain approximately constant. This 
effect supports the "downsizing" scenario of star formation in galaxies. While the general dependence 
of the galaxy color distributions on the environment is small, we find that the change of red galaxy 
fraction with epoch is a function of the local galaxy density, suggesting that the downsizing effect may 
operate with different timescales in regions of different galaxy densities. 
Subject headings: Galaxies: evolution — galaxies: fundamental parameters 



1. INTRODUCTION 

Many investigations in the past 30 years have shown 
that galaxy populations, manifested by galaxy colors, 
spectra, and morphological distributions, have a strong 
dependence on the environment and luminosity of the 
galaxies (e.g., Melnick & Sargent 1970; Dressier 1980; de 
Vaucouleurs 1961 and many other subsequent studies). 
More recently, much more detailed analyzes of galaxy 
color distributions in the local universe have become 
available using large galaxy samples of tens of thousands 
from the Sloan Digital Sky Survey (SDSS). Baldry et al. 
(2004) and others show that the colors of galaxies can be 
neatly separated into components of red and blue galax- 
ies. Balogh et al. (2004, hereafter B04) were able to 
delineate the dependence of galaxy colors on galaxy lu- 
minosity and environment, concluding that there is only 
a weak dependence on the latter. 

In this Letter we present a similar study of galaxy sam- 
ples at redshift between 0.2 and 0.6 using the photometric 
redshift galaxy catalogs from the Red-Sequence Cluster 
Survey (RCS) from Hsieh et al. (2005, hereafter H05). 
Combined with results from a Z'^0.05 sample from the 
SDSS (B04), we examine the trend of the dependence 
of galaxy colors on luminosity and environment at dif- 
ferent epochs. In §2 we briefly describe the photometric 
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redshift galaxy sample and the measurement techniques. 
We present the results in §3; their implications are dis- 
cussed and summarized in §4. We adopt a flat cosmology 
with = 0.3, = 0.7 and Hq = 70 km s'^Mpc'^ 

2. THE PHOTOMETRIC REDSHIFT GALAXY DATA 

The RCS is a 92 square degree imaging survey in the z' 
and Rc bands conducted with the CFHT 3.6m and the 
CTIO 4m to search for galaxy clusters in the redshift 
range of z < 1.4 (see Gladders & Yee 2005). Additional 
imaging in the V and B bands was also obtained for 
33.6 square degrees in the northern patches using the 
CFH12K camera. The observation and data reduction 
techniques are discussed in detail in Gladders & Yee 
(2005) for the z' and Rc data, and in H05 for the V 
and B data. Photometric redshifts for 1.2 million galax- 
ies using the four-color photometry are derived using an 
empirical training set method with 4,924 spectroscopic 
redshifts. Detailed descriptions of the photometric red- 
shift method and the precision and completeness of the 
sample are presented in H05. 

We select samples of galaxies covering two moderate 
redshift intervals using conservative criteria to minimize 
redshift and color errors. We choose galaxies with abso- 
lute magnitudes Mr^ < —18.5 in the photometric red- 
shift bins: 0.2 < z < 0.4 (hereafter, the z~0.3 sample) 
and 0.4 < z < 0.6 (hereafter, the ^'^0.5 sample) with 
redshift uncertainty fTz/i^ -I- z) < 0.1, where CTz is the 
computed photometric redshift uncertainty (see H05). 
The absolute magnitudes are corrected for K-correction 
and estimated evolution. For each galaxy, we use model 
colors computed from the spectral energy distributions 
of galaxies from Coleman et al. (1980) to estimate a 
galaxy spectral type based on the observed Rc — z' color, 
and derive K-corrections for the B and Rc magnitudes. 
We approximate the Rc band luminosity evolution by 
M{z) = M(0) - Qz (see Lin et al. 1999), where we 
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Fig. 1. — Upper large panel: Distributions of rest {B — R) colors for galaxies in subsamples of different surface density (S5) and luminosity 
i^Rc) 0-2 < 2 < 0.4. Subpanels for different densities are plotted horizontally, while subpanels for different absolute magnitudes are 
plotted vertically. Lower large panel: same for 0.4 < z < 0.6. Solid lines are two-Gaussian fits to the distributions of blue and red galaxies. 



the edges of the fields. These criteria produce a sam- 
ple of 106,095 and 124,004 galaxies for the z~0.3 and 
Z'^O.S sample, respectively. 

We estimate the local galaxy density using a projected 
surface density as a proxy (see, e.g.. Dressier 1980). For 
each galaxy in the primary sample, we measure i?5, the 
distance to the fifth nearest neighbor with M/j^ < —19.5, 
from which a local surface density parameter S5 in units 
of Mpc~^ (proportional to R^^) is computed. For count- 
ing the five objects we select only those with photometric 
redshifts within 2az of the primary galaxy. A complete- 
ness correction for each neighbor is applied, as not all 
galaxies within the sampling magnitude range satisfy the 
redshift uncertainty criterion. The completeness factor 
is estimated using the ratio of the total number of galax- 
ies within a, 0.1 mac bin at the ma.priitiide nf the neiprh- 



criterion. Since red galaxies on average have lower photo- 
metric redshift uncertainties, we refine our completeness 
correction by computing the factor separately for red and 
blue galaxies, separated at B — Rc — 1.8. The average 
completeness correction factor is ^^1.10. We count the 
nearest neighbors to the primary galaxy by summing the 
corrected weight for each galaxy until the total is > 5.0. 

To be able to put the E5 parameter for galaxies at 
different redshift on a comparative scale, a background 
correction must be applied to E5. This is estimated for 
each primary galaxy by computing the average number of 
galaxies in a circle of radius R5 using galaxies in the same 
RCS patch (of size ~ 2.5 deg 2, see Gladders & Yee 2005) 
that satisfy the redshift and magnitude criteria used for 
the counting of the neighbors. Completeness corrections 
are a.l.so annlied in the counting of the backaround ffalax- 
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Fig. 2. — Red galaxy fraction, fr, as a function of local galaxy 
density for samples of different luminosities for the two redshift 
bins. 
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Fig. 3. — Galaxy color {(B — Rc)o) peak of the blue and red 
galaxy distributions as a function of local galaxy density for sam- 
ples of different luminosities for the two redshift bins. 



of the background produces negative local densities; for 
convenience, we add 5 to S5 for the purpose of displaying 
some of the results. 

3. RESULTS 

The two panels of Figure 1 show the rest {B — Rc)o dis- 
tributions in bins of alDsolute magnitude and E5 for the 
2;~0.3 and bins. Wc divide the galaxies into sub- 

samples of local density such that the first bin contains 
the least dense 12.5% of the galaxies; the next three 
bins, each of the successive 25%; and the last bin, the 
most dense 12.5%. The solid lines are the two-Gaussian 
fits for the red and blue populations. Qualitatively, the 
two-Gaussian model that has been used for the SDSS 
z ~ 0.05 galaxy samples by Baldry et al. (2004) and B04 



and our data are: r ^ i?c+0.5 and u—r ~ 2(5 — i?c) — 0.4 
(based on galaxy color data from Fukugita et al. 1995). 
We also repeat the calculation using a sample with a 
more liberal error criterion, ^^/(l + z) < 0.4; this pro- 
duces qualitatively identical results, showing that the in- 
completeness correction works well. 

Figure 2 shows the fraction of red galaxies, fr, de- 
rived from integrating the Gaussian fits, as a function 
of local density for each magnitude bin for the two red- 
shift samples. The uncertainty for each fr value is esti- 
mated by generating 100 Monte-Carlo realizations based 
on the errors on the parameters of the corresponding 
two-Gaussian fit. We take the 68% confidence limits on 
fr for the 100 results to be the ilcr uncertainty. For 
both redshift bins, there is a strong dependence of fr on 
luminosity, covering a factor of 5-6, with the luminous 
galaxy samples having the higher values. Within each 
magnitude bin, the dependence of fr on local density, 
however, is relatively moderate, especially for the fainter 
galaxies. 

Figure 3 illustrates the dependence of the peak 

{B — Rc)o of the Gaussian fits on local density for galax- 
ies of different luminosities. Similar to the SDSS result 
at 2~0.05 (B04; Hogg et al. 2003), the peak colors of the 
color distributions of both the red and the blue galaxies 
have only a weak dependence on the local galaxy den- 
sity. We note that because of the scatter of the back- 
ground corrections, there will be some blurring of the 
local density parameter; however, this smoothing should 
have minimal effect at the high density bins, where the 
density is much higher than the scatter in the background 
correction. 

4. DISCUSSION AND CONCLUSIONS 

Using the three redshift epochs, we find a number of 
interesting evolutionary effects in the color distribution 

of galaxies. We find the peak colors of the red galaxy 
distributions to be remarkably similar over the three 
redshift bins. The color peaks for both the z^O.3 and 
samples for the different magnitude subsamples 
range from (B — Rc)o of ~ 1.3 to ~ 1.6, with the brighter 
red galaxies being redder. This is essentially identical to 
the z^0.05 SDSS samples of B04 over a similar luminos- 
ity range, which have u—r ^ 2.2 to 2.5 (or B—Rc ^ 1.3 to 
1.5). Hence, there has been only minimal changes in the 
colors of red galaxies in all environments from 2;~0.5 to 
z^0.05. The small color changes suggest that red galax- 
ies of all luminosities and in different environments are 
already well-evolved by redshift ^ 0.6 (e.g., see models 
in Bruzual & Chariot 1993). 

The difference in {B — Rc)o for red galaxies of different 
absolute magnitudes is likely due to the slope of the color- 
magnitude relation (CMR) for early-type galaxies. The 
range of (B — Rc)o of ^ 0.3 mag over ^ 4.5 mag in Mji^ 
is equivalent to a CMR slope of ~0.06, very similar to 
the slope of 0.03 to 0.08 found in the CMR of galaxy 
clusters (e.g., Lopez-Cruz et al. 2004). 

The blue galaxy color distributions, however, show a 
strong evolutionary effect from z^O.5 to z^O. For the 
z^O.5 sample, the peak colors of the blue galaxies are 
essentially independent of the luminosity of the galax- 
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blue as the faintest blue galaxies at lower redshift. For 
the 2:~0.3 sample, the blue galaxies in the two brighter 
bins {Mr^ < —20.5) now have color peaks 0.2 to 0.4 mag 
redder than that of the faintest blue galaxy bin. Com- 
paring to the 2~0.05 sample of B04 (their Figure 3), the 
evolution is quite dramatic. Here, the peak color of each 
successive brighter bin shifts significantly to the red, to 
the point where the brightest blue galaxy bin has almost 
the same color peak as the faintest red galaxy bin. We 
also note that the bluest peaks in our data are about 0.1 
mag bluer than the corresponding bin in B04. 

This evolution of the blue galaxy peak colors can be 
interpreted in the context of a 'down-sizing' scenario, a 
term first suggested by Cowic ct al. (1996). To first or- 
der, we can assume that (S — i?c)o provides a measure 
of the specific star formation rate (SSFR, defined as the 
star formation rate per imit mass) , or the time since the 
last major star formation episode (e.g., Bruzual & Char- 
lot 1993). At blue galaxies of all luminosities have 
very similar {B — Rc)o color peaks, which is consistent 
with galaxies of diff'erent masses having similar SSFRs. 
The mean colors of the brighter blue galaxies then be- 
come progressively redder with time, while those of the 
faint galaxies remain blue. This can be interpreted as 
the decrease in the SSFR being a strong function of the 
luminosity of the galaxy: the more massive the galaxy, 
the larger the decrease, or alternatively, the earlier in 
time the decrease takes place. 

This interpretation is in excellent agreement with 
the results of Heavens et al. (2004), who analysed the 
2;~0.05 SDSS spectra to recover the star formation his- 
tory. Their Figure 2 shows that low-mass galaxies have 
a more or less constant star formation rate from z ^ 1, 
peaking at z ^ 0.4 to 0.2. In comparison, the more mas- 
sive galaxies have their star formation peak much earlier, 
and show a decline as early as ^! ~ 1, and by .2~0.05 the 
SSFRs for galaxies of different masses are different by 
a factor of more than 10. Wc note the complementary 
nature of the techniques used to arrive at a similar con- 
clusion, in that the Heavens et al. work is based on the re- 
construction of star formation histories from present-day 
stellar populations, whereas our results are derived from 
direct measurements of galaxy colors at earlier epochs. 

We find clear indications of evolution in fr from 
to z^0.05, in that there is a general increase in 
fr at lower redshift. Comparing our high redshift bins 
with the SDSS results (B04), the evolution of fr also ap- 
pears to be a function of the local galaxy density: the 
2^^0.5 sample has a much flatter dependence of fr on 
local density. In B04's Figure 2 there is a strong local 
density dependence of fr for galaxies of all luminosities; 
fr ranges from 0.1 to 0.4 for galaxies of different lumi- 
nosities in the lowest density sample, and from 0.65 to 
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0.85 for the highest density samples. By comparison, fr 
in our two higher redshift bins has a similar range (0.1 
to 0.4) for the low density environment, but has values 
between 0.15 and 0.6 in the densest regions. Thus, there 
is a lack of faint red galaxies in the high density bins: 
fr < 0.4 for all but the most luminous galaxies. In other 
words, at z^O.2, there arc many moderately luminous to 
faint blue galaxies in high density regions. This drop of 
fr at high density and high redshift could be a reflec- 
tion of the change in the morphology-density relation at 
z ~ 0.5 (e.g.. Dressier et al. 1997; Smith et al. 2005). 
The lack of faint red galaxies in high density regions has 
also been noted by others at moderate to high redshifts 
(e.g., Kodama et al. 2004) 

The density-dependent evolution of fr suggests the 
need for an additional element in the downsizing sce- 
nario: the downsizing timescale is a function of the en- 
vironment. That is, for galaxies of similar masses, the 
drop in star formation rate occurs successively later for 
galaxies in less dense environments. The increase of the 
faint blue galaxy fraction in the high redshift, high den- 
sity samples can be interpreted as these galaxies having 
their peak star formation era at a redshift of ~ 0.3 to 
0.2 (Heavens et al. 2004), followed by a decline, resulting 
in a high faint red galaxy fraction at 2:~0.05 (Figure 2, 
B04). In the less dense regions, the star formation rate 
continues to be high to lower redshifts (see, e.g., Gomez 
et al. 2003), keeping the blue galaxy fraction high. This 
is supported qualitatively by the positive, if small, gra- 
dients in the {B — i?c)o peak going from the least to the 
most dense samples, for all redshift bins. 

In summary we have analysed the {B — Rc)o color dis- 
tribution of a sample of ^ 230, 000 galaxies with photo- 
metric redshifts between 0.2 and 0.6 and Mjj^ < —18.5, 
in combination with the results derived at z~0.05 by B04 
using SDSS data. We find there is a strong evolution in 
the peak {B — Rc)o colors of the blue galaxies, consistent 
with a downsizing star formation scenario where more 
massive galaxies have their peak star formation time oc- 
curring at earlier epochs. This is in excellent qualitative 
agreement with star formation histories of galaxies of dif- 
ferent masses derived by Heavens et al. (2004) using the 
technique of "fossil stellar populations." We also find 
that the evolution of the red galaxy fraction is depen- 
dent on local galaxy density, in that there is a significant 
deficit of faint red galaxies in rich environments at high 
redshift. In a future paper we will examine the evolu- 
tion of the color distributions for galaxies in different 
local densities in various global environments (e.g., clus- 
ter core, cluster outskirt, and field), and interpret the 
results more quantitatively in terms of stellar mass and 
star formation rate via stellar population models. 
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